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ABSTRACT 


For a plane wave incident on a perfectly conducting ephere of radius a, 
the surface field components are expressible in terms of series of the Mie type. 
These have been evaluated on an IBM 7090 computer and results are presented 
for @=0(5)180°, where 6 denotes the angular position on the surface, and 


ka = 0, 1(0, 1)5. (0, 2)10.0, where k is the propagation constant, 
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INTRODUCTION 

The increasing interest in sophisticated methods of controlling the scatter- 
ing behavior of bodies using, for example, absorbers and reactive loading, has 
served to focus attention on the surface field, a knowledge of which is almost 
essential if these techniques are to be employed in anything near an optimum 
manner. Although the ultimate requirement is still, in general, the determina- 
tion of the far field, it is now appreciated that an examination of the surface field 
can provide new insights into the scattering processes, An assumption about the 
currents induced on the surface of the body ix often the basis for approximate tech- 
niques for cross section estimation, and the surface field is also the natural pro- 
duct of digital methods of far field computation involving the numerical solution of 
the integral equations (see, for example, Andreasen, 1965) or the point-by-point 
satisfaction of the boundary conditions (Mullin et al, 1965). 

Studies of the surface field either by probe measurements (Senior, 1964) or 
by digital techniques have made it desirable to have some body whose field is known 
precisely and which can be used for calibration or checking purposes, The sphere 
is the natural candidate for this role, For plane wave incidence, an exact expres- 
sion for the scattered field at all points of space is available in the form of an 
infinite series involving spherical Bessel and Legendre functions, The series has 
been quite extensively computed in the far zone, particularly for the back scatter- 


ing direction (Hey et al, 1956; Bechtel, 1962), but curves showing the bistatic 
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scattering in the E- and H-planes ave also been given by King and Wu (1959) for 
17 values of ka in the range 1,1 to 20, By comparison, however, the expressions 
for the surface field have received less attention, and it would appear that only 
plots for isolated values of ka have been obtained, 

In the course of a recent investigation of reactive loading :.pplied to a 
metallic sphere (Liepa and Senior, 1964), it was necessary to calculate the surface 
field components for a wide variety of ka and 6, where @ is the angular position 
on the sphere, The series were programmed for the IBM 7090 computer of The 
University of Michigan, and it seemed worthwhile to carry through the computation 
for a reasonably comprehensive set of values to permit a detailed numerical and 
conceptual picture of the surface field behavior, The tables presented here give 


both tangential components of the magnetic vector for ka=0. 1(0, 1)5, 0(0, 2)10, 0 


and @=0(5)180°, and required 10,5 minutes of machine time to compile, 
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EXPRESSIONS FOR THE SURFACE FIELD COMPONENTS 
Consider a perfectly conducting sphere of radius a whose center is at the 
origin of a Cartesian coordinate system (x,y,z). A plane electromagnetic wave 
is assumed incident in the direction of the negative z axis and since there is no 
loss of generality in taking its electric vector to lie in the x direction, we choose 
i 


Eb ote . H = -yYe (1) 


where Y is the intrinsic admittance of free space and a time factor eo has been 
suppressed, 
if we also introduce the spherical] polar coordinates (r,@,$) such that 


x=rsinOcosf, y=rsinO sing, z=rcosé@ 


with @=0 representing the back scattering direction and @=7 the forward one, 
expressions for the total (incident plus scattered) field can be obtained from the 
standard Mie solution. For the magnetic vector in particular we have, on the 
surface r=a, 


H_=0 
r 


Be = 6 n+i2n+1 aati {a 1 Pare 4 i i, «| K r) 
7 = sin Di n(n) ion (ka) = sin@ ray 7) oe 
(2) 


Zi i pt (cos 8) 
nicos ont g ka) sin@ 


(3) 
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(see, for example, Liepa and Senior, 1964), where 


¢ (ka) = ka (ka) 


and ba) is the spherical Hankel function of the second kind, The primes 
denote differentiation with respect to ka, P)(cos 6) is the Legendre function of 


degree n and order unity as defined, for example, by Stratton (1941). 


The tangential components H, and Hy are directly related to the surface 


6 


current J via the equation 
J=%xH , 
implying 
J,=-H,, J,=H, , 
6 f°’ ‘$ 


and for convenience they are written as 


Ho = Ysing T, (6) 2 


Hy = Y sing T,(6) . 


T, and T, are, of course, functions of ka as wellas 9, but are independent of 


2 


~. From equations (2) and (3), their expressions are 


1 
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T(@) = is = + atari) 


ti(ka) sind © (ka) 00 oP Soma (4) 


1 
: P (cos 6) 
ae ntl 2ntl 1 i n 
T(6) ” ka DB : ao a ¢'(ka) 7 3” “(cos 6) * © (ka) Cn t (ka) ~ sin6 } (5) 


and these are the functions that were computed, 
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COMPUTATION OF T, (6) AND T,(6) 


The series representations ot T and Ty shown in equations (4) and (5) 
were programmed for The University of Michigan IBM 7090 computer, The pro- 
cedure was quite straightforward, Because of the similarity of T, and Ty the 
terms in both series were evaluated within the same loop, and the series were 


terminated whenever the magnitudes of successive terms fell below 107". With 
this criterion, the number of terms included was, in general, the same for Tv 


and T,. It depended primarily on the value of ka and increased from 4 for 
ka=0.1, through 9 and 17 for ka=1 and 5 respectively, to 25 for ka=10, 

The evaluation of the spherical Hankel and Legendre functions was carried 
out using external subroutines which were originally prepared in connection with 
an investigation of reactive loading applied to a sphere, Although the present com- 
putation was limited to ka < 10, and required only the first 25 (cr fewer) terms in 
the series, the subroutines are sufficient to produce at least five digit accuracy 
for n upto 44, The program is therefore capable of cciaputing T, and Ty for 
ka as large as (say) 20 with no appreciable modification. 

A few comments about the two subroutines may be of interest, The function 


7 (ica) was written as a sum of spherical Hankel functions in the form 


- 1 - 
e"(x) = arf coe, + h ix , 


and h (x) was itself broken up into spherical Bessel functions of the first and 
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second kinds: 


19 « jt Od 


n 


j,( was obtained by numerical integration of the finite integral expression 


G) 


n, 


j,09 = 


n[2 
\ cos (xsin pcos z dg 
0 


(see, for example, Adams and Hippisley, 1947). The range of integration was sub- 
divided into 40+2n intervals, and to iudge from spot checks, the resulting evalu- 
ations were accurate to six significant figures for n< 20, and to five for 20<n< 44, 


The spherical Bessel functions of the second kind were determined from the finite 


series exparsion 


nt+1 r ' 
y (x) = cy {= (x+7) Stee or 
r=0 (2r)! (n-2r)! (2x) 


: nt (-1)* (n+2r+1)! 
= sin (x+ a) ue ae Cee 
r=0 (2r+1)! (n-2r-1)! (2x) 


(Watson, 1948), giving seven digit accuracy for n < 20, and five or better to n=44, 


The most direct method of computing the Legendre function factors is to use 


the recurrence relations 


1 _ 2nt+1 1 nti i 
Pp nv (008 9) a cos@ P (cos 4) os P1008 9) 
1 1 
P (cos@) P” .(cosé@) 
~2_ p!(coso) = ncos§ — ~ (n+1) al 
36 on sin sin§ 
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but because of the loss of accuracy accompanying their repeated use, equations (6) 
and (7) were employed only for the generation of functions for which n >32. For 
n< 32, Tallgvist (1938) bas given the expansion of P (cos6) in terms of-cosines of 


multiple angles, and since 
P*(cos6) = ©. P (cose) 
n ~ 80°n : 


direct differentiation leads to a comparable «<pansion of P*(cos6) in terms of sines 
of multiple angles. This was used to compute P’(cos6) for n< 32, with 

# P*(cos6) determined from equation (7). The values of the two Legendre functio 
obtained in this manner are believed ccurate to seven significant figures for 


O0<¢n<¢32, and to five significant figures for 32 <n < 50. 


THE UNIVERSITY OF MICHIGAN 
59548-3-T 


GRAPHS AND TABLES 

The functions T,() and T,(6) were computed for ka=0, 1(G, 1)5. 0(0. 2)10. 0 
and for 0< @ < 180° in increments of 5°, and in the following tables the real and 
imaginary parts are presented along with the magnitudes and arguments, The 
arguments are given in degrees and since these are, of course, principal values, 
the addition (or substraction) of integer multiples of 360° may be necessary to 
convert them into physically-meaningful phase angles, In this connection we note 
that the incident field was chosen to have zero phase in the plane of the shadow 
boundary. We also remark that for ka=0, T, (6) =- . cos@ and T,(6) =~ 5 ; 

No detailed study of the accuracy of the tabulated values has been carried 
out, but from spot checks using a hand computer and standard tabulations of the 
Hankel and Legendre functions, it is believed that the results are accurate to all 
of the decimal places show . 

To illustrate the nature of the surface field components, the amplitudes 
and phases of T,@) and T,(6) are plotted as functions of @ for ka=0,1, 1,0, 
5.0 and 10.0 in Figures 1 through 4 respectively. The values were taken directly 
from the tables except for the addition of appropriate multiples of 360° to the 


listed arguments in accordance with the known physical behavior of the field, 
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FIG, 1; SURFACE CURRENT COMPONENTS FOR ka=0,1 
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FIG. 2; SURFACE CURRENT COMPONENTS FOR ka = 1.0 
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FIG, 3; SURFACE CURRENT COMPONENTS FOR ka=5,0 
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FIG, 4. SURFACE CURRENT COMPONENTS FOR ka= 10," 
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